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Abstract 

It is shown by numerical calculations that the convoluted QCD pomeron propagator in the 
external field created by a solution of the Balitsky-Kovchegov equation in the nuclear matter 
vanishes at high rapidities. This may open a possibility to apply the perturbative approach for the 
calculation of pomeron loops. 

1 Introduction 

In the QCD, in the limit of large number of colours, strong interaction at high energies is mediated 
by the exchange of BFKL pomerons, which interact via their splitting and fusion. In the quasi- 
classical approximation for photon (hadron)-nucleus scattering the relevant tree (fan) diagrams are 
summed by the well-known Balitsky-Kovchegov (BK) evolution equation [1, 2, 3]. For nucleus- 
nucleus scattering appropriate quasi-classical equations were derived in [4, 5]. In both cases 
pomeron loops were neglected. This approximation can be justified if the parameter 7 = A exp Apy 
is small, with y the rapidity and Ap and A the pomeron intercept and triple pomeron coupling. 
Then for a large nuclear target, such that A^^'^-y ~ 1, the tree diagrams indeed give the dominant 
contribution and loops can be dropped. However with the growth of y the loop contribution becomes 
not small and this approximation breaks down. 

Direct calculation of the loop contribution seems to be a formidable task for the non-local 
BFKL pomeron. Simplest loops have been studied in several papers for purely hadronic scattering 
[6, 7, 8]. In particular in [8] it has been found that pomeron loops become essential already at 
rapidities of the order 10-^15. They shift the position of the pomeron pole to the complex plane and 
thus lead to oscillations in cross-sections. However with the growth of energy loop contributions 
begin to dominate and one needs to sum all of them. There have been many attempts to do this 
in the framework of the so-called reaction-diffusion formulation of the QCD dynamics and the 
following correspondence with the statistical approach [9, 10, 11, 12, 13, 14] (see also a review 
[15] and references therein). Unfortunately concrete results could be obtained only with very 
crude approximations for the basic BFKL interaction and the stochastical noise in the statistical 
formulation. The conclusions of different groups are incomplete and contradictory. So in [14] it was 
found that the geometric scaling following from the BK equation was preserved with loops taken 
into account, although going to the black disc limit was much slower. On the contrary in papers 
based on the analogy with statistical physics (see [13, 15]) it was argued that the BK scaling was 
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changed to the so called diffusive scaling (with an extra in the denominator of the argument) 
but the speed of achieving the black disk limit was essentialh'' unchanged. 

In our previous study of pomeron loops [16] we considered a much simpler model with the 
local supercritical pomeron in the Regge-Gribov formalism. Instead of trying to solve the model 
for the purely hadronic scattering we considered the hadron-nucleus scattering and propagation 
of the pomeron inside the heavy nucleus target. Moreover to avoid using numerical solution of 
the tree diagrams contribution with diffusion in the impact parameter, we concentrated on the 
case of a constant nuclear density, which allowed to start with the known analytical solutions. We 
have found that the nuclear surrounding transforms the pomeron from the supercritical one with 
intercept e > to a subcritical one with the intercept — e. Then Regge cuts, corresponding to loop 
diagrams, start at branch points located to the left of the pomeron pole and their contribution is 
subdominant at high energies. As a result, the theory acquired the properties similar to the Regge- 
Gribov one with a subcritical pomeron and allows for application of the perturbation theory. In [16] 
we expressed our hopes that a similar phenomenon might occur in the QCD with BFKL pomerons. 

In this note we demonstrate that such hopes are possibly founded. We consider the pomeron 
propagator in the nuclear field and give arguments that, similar to the local Regge-Gribov case, 
it vanishes at large rapidity distances. We stress that at present we are unable to give the full 
proof for this behaviour. Our study is based on numerical calculations. This makes us to choose a 
relatively small subset of initial conditions out of the complete set necessary for the study of the 
pomeron propagator. Our numerical results show that, with the chosen set of initial conditions, 
this pomeron propagator vanishes at large rapidity distances. This result may serve as a starting 
point for the study of loops and can be directly applied for double inclusive cross-sections in the 
nucleus-nucleus scattering. 

2 Local pomeron in the nuclear field 

To make clear our motivation and goals we start with reviewing our treatment of a simpler model, 
the local reggeon field theory (LRFT) with a supercritical pomeron and triple pomeron interac- 
tion. The following passage to the non-local BFKL pomeron will not introduce ideologically new 
moments, although greatly complicates the technique. 

The LRFT model is based on two pomeron fields b) and b) depending on the rapidity 
y and impact parameter b, with a Lagrangian density [17] 

L = Lo + A$^$($ + + gp^K (1) 

Here the free Lagrangian density is 

Lo = dy -a'Vl + e)$ = ^^S^, (2) 

where e is the intercept minus unity and a' is the slope. The source term describing interaction 
with the nuclear target at low energies is 

gp{y,b)=gATA{b)6{y), (3) 

where g is the pomeron-nucleon coupling constant and T{b) the profile function of the nucleus. For 
a supercritical pomeron e > and A < 0. 
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Figure 1: The new vertex for twopomeron annihilation, which appears after the shift in field $ 



Solution of the classical equation of motion for field is $o = and the classical equation for 
$ becomes 

dy^ = -a'Vl^ + e$ + A$^ (4) 

with an initial condition 

^y = 0)=gTA{b). (5) 

Equation (4) describes evolution of the pomcron field in rapidity and its diffusion in the impact 
parameter inside the nucleus. We denote the solution of the classical equation of motion (4) with 
the initial condition (5) as ^o{y,b). 

To go beyond the classical approximation and thus study loops one makes a shift in the quantum 
field 

^y,b) = My,b) + ^i{y,b) (6) 

and reinterprets the theory in terms of quantum fields $i and In the Lagrangian terms linear 
in vanish due to the equation of motion for $ and we obtain 

L = + 2A$o)^^ + A$o^^^ + A$i$t($i + $^). (7) 

This Lagrangian corresponds to a theory in the vacuum with the pomeron propagator in the external 
field /(6,y) = 2A$o(2/,&) 

P = -(S + 2A$o)~\ (8) 

the standard triple interaction and an extra interaction described by the term A$o^^^- This new 
interaction corresponds to transition of a pair of pomerons into the vacuum at point (y, b) with 
a vertex X^Q{y,b), see Fig. 1. The propagator P in the external field is described by a sum of 
diagrams shown in Fig. 2. 

Loops can be formed both by the standard interaction and the new one. In the latter case they 
are to be accompanied by at least a pair of standard interactions. Simplest loops in the Green 
function are illustrated in Fig. 3, where it is assumed that the propagators are all taken in the 
nuclear field /. Note that a loop formed by the standard interaction has the lowest order A^/a' 
and requires renormalization in this model. A loop formed by the new interaction has the lowest 
order X^/a' and is finite. 

The quantum part of the amplitude is obtained as a tadpole g < $i(y, 6) >. The simplest 
diagram for it contains one loop and is shown in Fig. 4. 

Propagator (8) in the external field f{y,b) = 2\^o{y,b) satisfies the equation 

^^^y^^^^y'^'^ ={e- a'Vl)P{y, b\y', b') + f{y, b)P{y, b\y' , b') (9) 
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Figure 2: Diagrams summed into the pomeron propagator in the nuclear field 
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Figure 3: Simplest loop diagrams for the pomeron Green function 




Figure 4: The lowest order diagram for the quantum correction to the scattering amplitude 
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with the boundary conditions 



P(2/,6|y',60 = O, y-y'<0, P{y' ,b\y' ,b') = 5^{b - b'). (10) 

In the general case propagator P can only be calculated numerically, just as the external field 
f{y, b). Its analytic form can be found in two cases [16]: if the slope a' = or if one is studying 
the nuclear matter with the profile function TA{b) = Tq independent of h. In the second case, 
relevant for our following studies, the nuclear field becomes independent of h (but dependent on y) 
and the propagator in this field takes the form 

1 P^{y') Ay-y')- 



P{y,b\y,b') = —^, L^e^yy-y-T^r^^)^ (n) 

where 

P{y) = l-^{ey^-l). (12) 
Remarkably at y ^ oo this propagator vanishes as e' 

-{y-y')e_ This means that contributions from 
the loops will be small at large y, so that one can calculate them by the standard perturbation 
approach. 

3 Interacting BFKL pomerons 

Now wc pass to our main subject, the theory of interacting BFKL pomerons. By its structure it 
is quite similar to the considered LRFT model. However now the pomeron, apart of the impact 
parameter b, depends on the relative transversal distance r between the two reggeized gluons out 
of which it is formed. So it is described by two fields ^{y,ri,r2) and ^*^(y,ri,r2) depending on 
rapidity y and transverse coordinates of the reggeized gluons ri and r2- The Lagrangian has a 
general structure similar to (1): 

L = Lo + Li + Le. (13) 

Here the free Lagrangian is 

/■ f) 
Lo = J dVidV2$^ V? Vi (— + Hbfkl)^, (14) 

where Hbfkl is the Hamiltonian for the BFKL pomeron [18]. The triple pomeron interaction is 
given by 



Lj = ^^ I ''I l\'^^ ^\y,n,r2)^\y,r2,T^)K^My.r^,n) + o $t), (15) 
where rik = r^ — and 

K^i = r-liViVi (16) 
The interaction with the nuclear target is 

Le = - j d^riSr2^U, (17) 

where 

J{y,rur2) = gATA{b)p{ri2)6{y), (18) 
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Ta is as before the profile function of the nucleus and p is the colour density of the nucleon. 

Our idea is to repeat the procedure applied to the local pomeron as described in the preceding 
section. First one solves the classical equation for the fields. Its solution for $ is a sum of all tree 
diagrams with the pomeron splitting in two. This solution satisfies the BK equation and although 
not known analytically can be found numerically. At the next step we calculate the quantum 
correction to this solution (pomeron loops). To this aim wc shift the field $ with a non-zero 
classical value $o by this value, passing to the new quantum field $i = $ — exactly as for the 
LRFT, and study the theory in terms of fields $i and This leads to a change in the Lagrangian 
quite similar to what has been described in the preceding section. In the new theory the propagator 
is to be taken in the nuclear field. This is achieved by substituting one of the two fields $ by its 
classical value $o in tlie part of interaction quadratic in $. Also in the new theory there appears 
a new interaction term, obtained when in the part of interaction linear in $ one substitutes $ by 
its classical value. 

In this paper we are not attempting at the calculation of loops. We limit ourselves with the 
starting point of this program: studying the pomeron propagator in the described nuclear field. 
Our aim is to understand if this propagator vanishes at large rapidities instead of growing in the 
absence of the field, as it occurs in the LRFT. If it is indeed so then one can hope to calculate 
loop contributions in the nuclear surrounding perturbatively. Otherwise this surrounding is useless 
for loop calculations. The BFKL propagator P{y,ri,r2;y' ,r'i,r'2) in the nuclear field satisfies the 
equation 



dP{y,ri,r2) " /" ,2 ^ 
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dy 27r J rfgrjg 

{P{y,ri,n) + P{y,r2,rs) - P{y,ri,r2) - ^{y,ri,rz)P{y,r2,r^) - $(y, r2, r3)P(y, ri, ra)) , (19) 

where a = agNd'K and we suppressed the second three variables coordinates y', , r2 on which the 
equation does not depend. The dependence on them follows from the initial condition: 

P{y = y',ri,r2;?/',ri,r^) = V^''V^H\n - r[)5\r2 - r'^). (20) 

Function $(y, ri,r2) in this equation is the solution of the BK equation for the sum of all fan 
diagrams in the nucleus. It satisfies a non-linear equation 

— 1 =7r d'^r3^^[^{y,ri,r3) + ^{y,r2,r3) - ^{y,ri,r2) - ^iy,ri,r3)^{y,r2,r3)) 

Oy ZTT J ^13''23 

(21) 

with the initial condition appropriate for a given nucleus 

$(2/ = 0,ri,r2) = $o(ri,r2) = 1 - e-^^^^^P^''^^\ (22) 

We remind that h= {r\ + r2)/2, Ta is the conventional profile function of the nucleus and p is the 
colour density of the nucleon. 

We shall again study the simple case of the nuclear matter when Tyi(6) = Tq is independent of 
h. Then from the structure of Eq. (21) it follows that $ is a function of only the relative distance 
between the gluons ^•(j/, ri,r2) = $(y,ri2). Since the study is only possible numerically, to avoid 
using the singular initial condition, we shall consider a convolution of P(y, ri,r2; y', , ^2) with an 
arbitrary function Vi'V2ip{ri,r2)'- 

P{y,ri,r2) = I d\[d'r'2Piy,n,r2-,y',r[,r'2)VMijir[,r'2). (23) 
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This convolution satisfies the same equation as the propagator itself but at y = y' we have 

P{y = y',ri,r2)=4^{n,r2). (24) 

Obviously the properties of the propagator can be studied taking a full set of functions V'(^i) ^^2)- 

Note that in the following, for simplicity, we mostly set y' = 0. The dependence on the choice 
of y' will be studied for the forward case and shown to preserve our main conclusions. 



4 Forward case 
4.1 Main equations 

Calculation become especially simple for the forward case when not only $ but also the convoluted 
propagator P depend only on the relative distance between the gluons in the pomeron. In this case 
the equation for P simplifies to 

dP{y, x) ^ f i2 ^'^ 



dy 



^ f (fxi^(p{y,xi) + P{y,X2)-P{y,x)-2^y,xi)P{y,X2)). (25) 

Here we denoted x = ri2, xi = ri3 and X2 = with X1 + X2 = x. In the same notation the nuclear 
field $ satisfies 



d^iy,x) 
dy 



- p 2 

= 7r / d^xi^^(Hy,xi) + ^iy,X2) - Hy,x) - ^y,xi)^y,X2)), (26) 

ZTT J XIX2 ^ ' 

We have to solve this pair of equations with the initial condition for $(?/ = 0, x) fixed by the 
properties of the nuclear medium, and the initial condition (24) for P{y,x) taken arbitrary, since 
we are interested in the propagator in a given nuclear surrounding. 

Equation (25) is a linear equation for P{y, x) in contrast to the BK equation. At y — t- 00 
^{y,x) — >■ 1 independent of the chosen initial condition. One may think that at y ^ 00 the 
behaviour of P{y, x) can be derived from the asymptotic equation 

^^^^'""^ =^ f (fx,^(p{y,x,) + Piy,X2)-Piy,x)-2P{y,X2)) 



dy 



= -Ffe,x)|/dV;||. (27) 

However the integral on the right-hand side has become divergent (although it converges at finite 
y). This means that the limit y — > 00 is more delicate and cannot be taken under the sign of integral 
over xi. And indeed we shall see by numerical calculations that the behaviour of the solution at 
y — >■ 00 is not solely determined by the limiting value of ^{y,x) but depends on its behaviour at 
finite y. 

For numerical studies both the BK equation and linear equation (25) in the momentum space 
are more convenient. Introducing 

,/ N '^iy,x) , . P{y,x) 

ny,x) = — Piy^^) = — 2 — (2^) 

and then passing to the momentum space we obtain the following equations for (p{y, k) and p{y, k) 

= -a{HBFKLHy, k) + 4>\y, k)) (29) 
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Figure 5: Solution (p{y, q) of Eq. (29) as a function of for different of y = ay. Curves from 
bottom to top correspond to y = 0, 2, 4, 6, 8 and 10. 



and 



where 



= -a(HBFKL + 2</.(y, k))p{y, k), 



Hbfkl = lnk'^ + lnx^ - 2(VXl) + ln2). 



(30) 



(31) 



To study the behaviour of the propagator in the external field (p one has to solve this pair of 
equations with the initial conditions 

<P{y, k)y=o = Mk), P{y, k)y=o = po{k), (32) 
with some fixed (po and for a complete set of function po{k). 

4.2 Numerical studies. Momentum space 

Wc have set up a program which simultaneously solves the pair of equations (29) and (30) for a 
given pair of initial conditions (32) by the Runge-Kutta method. For the BK evolution we have 
fixed the initial condition as 

(f)Q[k) = -(l/2)Ei(-A;V0.3657) (33) 

used in our previous calculations. We have limited our momenta in the interval 10~^ < k < 10^ 
and divided this interval into 1600 points. 

The behaviour of 0(A;) with k'^ at different values of the scaled rapidity F = ay = 2, 4, 6, 8 and 
10 is shown in Fig. 5. (Note that the maximal value of the scaled rapidity y = 10 corresponds to 
the natural rapidity of order 50). 

For the BFKL evolution in the external field </>, in the first run (A), we have taken the same 
form of the initial condition at y = but with a variable slope 



Po{k) = -(l/2)Ei(-A:Va). 



(34) 
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Figure 6: Solutions p{y, q) of Eq. (30) as a function of q'^ for different Y = ay for run A with a = 1. 
Curves which start to faU at higher q^ correspond to higher Y = 0,2, 4, 6, 8 and 10 

We have performed calculations for a = 0.2, 0.6, 1.0, 1.4 and 1.8. In the second run (B) the initial 
condition was taken with extra powers of k^. 

po(fc) = -(l/2)fc2"Ei(-fcV0.3657) (35) 

with n = 0, 1, 2, 3 and 4. 

In all cases the behaviour of the solution p{y, k) was found to be universal. At large enough y 
the solution becomes independent of k"^ up to a certain maximal ^^^^.(y), starting from which it 
goes to zero. Roughly 

p{y,k)^A{y)e{kl,,-k'). (36) 

As y grows A{y) goes to zero and k^riaxiv) go^s to infinity. So on the whole the solution vanishes 
as y — )• oo, its X dependence tending to 6'^{x). 

We illustrate this behaviour in Figs. 6 and 7, in which we show the solution p{y, k) for run A 
with a = 1.0 and run B with n = 2 as a function of k"^. One observes that although the values of 
p{y, k) for the two cases are different, their behavior with y is the same: they vanish as y — t- oo. 

This universality is especially obvious if one calculates the slope A(y, k) of the y-dependence of 
p{y, k) at fixed k presenting 

p{y,k)o^e^^^y^^\ (37) 

It turns out that at y > 1 the slope A(y,A;) is independent of k and identical for all considered 
cases (run A with all studied a and run B with all studied n). Its smooth behaviour with Y is 
shown in Fig. 8. One observes that starting from Y = h the slope becomes negative indicating that 
the solution goes to zero at y >> 1. 

One has to take into account that in the external field (j){y, k) depending on rapidity the pomeron 
propagator ceases to depend only on the rapidity difference. Rather the initial and final rapidities 
become two independent variables. To see what influence it has on the behaviour of the propagator 
at large rapidities we varied the initial rapidity y = y' for the evolution of p{y, k), leaving unchanged 
the initial rapidity y = for the evolution of ^(y, k), which is the rapidity of the nucleus. One finds 
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Figure 7: Solutions p{y, q) of Eq. (30) as a function of for different Y = ay for run B with n = 2. 
Curves which start to fall at higher q^ correspond to higher Y = 0,2, 4, 6, 8 and 10 
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Figure 8: Slope A of the F-dependence of the solutions p{y,q) of Eq. (30). The slope is the same 
for all initial conditions (34) and (35). 
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Figure 9: Solutions p{y, q) of Eq. (30) as a function of q'^ for different Y = ay for the initial value 
taken at Y = 3 according to run A with a = 1. Curves which start to fall at higher q^ correspond 
to higher y = 4, 6, 8 and 10 



that although at initial stages of evolution the behavior of p{y, k) strongly depends on the value 
of y', at higher rapidities this behaviour is essentially the same for any y', namely the convoluted 
propagator goes down with rapidity with the slope independent of yo- This is illustrated in Figs. 
9 and 10 which show results for y' = 3/a {Y' = 3) In Fig. 9 we show the solution p{y, k) for run 
A with a = 1. One observes, that although absolute values of p{y, k) are quite different from the 
case y' = shown in Fig 6, the behaviour with the growth of rapidity is the same. It is especially 
clear from the values for the slope A shown in Fig. 10 together with those for the case y = (Fig. 
8). Again at the initial stage of the evolution the behavior with y' = 3 is quite different from that 
with Yq = 0. However at higher rapidities the values for the slope become the same. 

It is remarkable that this behaviour takes place only with cj) given by the exact solution of the 
BK equation. Taking an approximate form 

$(y,x) ~ 1 -e-^'(2/)^', (38) 

where the "saturation momentum" Q^(y) ~ ^'^-O^ay -^q obtain an equation for p{y,k) in the mo- 
mentum space 

dp{y,k) _r^, „./ k^ 



-a 



HBFKL-Ei{-J^)]p{y,k). (39) 



dy L V 4Q2(y), 

Taking for simplicity the initial condition po{k) = 4>oik) at y = we get the solution shown in Fig. 
11. One observes that at large y the solution acquires the same form (36) where however A{y) 
grows with y: 

Aiy) ~ e^-^^. (40) 

This implies that with the approximate form (38) of $ the final solution P{y, x) in the coordinate 
space behaves in a singular manner at y — >■ oo. Effectively 

Piy,x)y^oo^e'-^yxH\x) (41) 

and it is impossible to say that it vanishes in this limit. 
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Figure 10: Slope A of the y-dependcncc of the solutions p{y, q) of Eq. (30) with the initial condition 
put at y = (upper curve) and at y = 3 (lower curve) according to run A with a = 1. 
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Figure 11: Solutions of Eq. (39) as a function of k'^ for different y. Curves from bottom to top 
correspond to Y = ay = 2, 4, 6, 8 and 10. 
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Figure 12: Solutions of Eq. (25) in the nuclear field as a function of for different y. Curves with 
a maximum from right to left correspond to Y = ay = 0, 2, 4, 6, 8 and 10. 



4.3 Numerical studies. Coordinate space 

It is instructive to compare the behaviour of the propagator in the momentum space with that 
in the coordinate space. Rather than perform the Fourier transformation of the found propagator 
p{y, k) and determine P{y, x) according to (28) we directly solved the forward equations in the 
coordinate space (25) and (26) choosing the initial conditions in accordance with (33) and (34). 

$(y = 0,a;) = l-e-3°=^^5^^' P(y = 0, x) = 1 - e-^"^'. (42) 

Again we used the Runge-Kutta method with 10~* < r < 10*. The number of divisions was taken 
400 in Inr^ and 200 in the azimuthal angle. 

Numerical calculations become much more cumbersome in the coordinate space due to the fact 
that the non-linear term and the term with the nuclear field depend on two different points. Having 
in mind that these calculations serve only for the illustrative purpose, we have restricted ourselves 
to only one value a = 0.3657 in (42), that is have taken the initial values for P{y, x) at y' = the 
same as for ^{y,x). Our results are shown in Fig. 12. The form of the solution in the coordinate 
space very much resembles that of the gluon density in the momentum space: there is a sharp 
maximum at a certain value of which with the growth of y moves in the direction of smaller 

(to the ultraviolet). However in contrast to the gluon density the height of the maximum is 
diminishing with y, so that in the limit y — >■ oo the propagator vanishes. 

In the construction of the amplitudes pomeron propagators are convoluted with the triple 
pomeron vertex. The behaviour of this convolution with the growth of rapidities can be studied 
from the integral 

hiy) = [ ^^^Pi{y,xi)P2{y,x2)P3{y,xi+x2) (43) 
where Pi,2,3 may correspond to the same or different initial conditions. Presenting Ii as 

/i(y) = aie-^i(f)^, y = aY (44) 
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Figure 13: Exponent Ai in (44) as a function of rapidity for the convolution of solutions to Eq. 
(25) with the initial conditions in the form (42). The curves from top to bottom correspond to 
1) all three solutions with a = 0.3657, 2) two with a = 0.3657 and one with a = 1.4, 3) two with 
a = 1.4 and one with a = 0.3657 and the lowermost to all three with a = 1.4. 

we find values of the exponent Ai from our numerical calculations shown in Fig. 13 for the initial 
conditions (42) with a = 0.3657 and a = 1.4. As one observes, although at lower rapidities Y the 
values of Ai depend on the choice of the initial conditions, starting from Y = 6 this dependence 
disappears and the behaviour of Ai becomes universal. In any case it is clear that the convolution 
vanishes at F — > oo. To compare, it grows as exp{3ABFKLY) in the vacuum, with Abfkl = 41n2. 

5 General case 

Calculations of the loop contribution require the propagator in the non-forward direction, depending 
not only on the relative distance between the gluons but also on the center-of mass coordinate. The 
convoluted propagator P{y,ri,r2) in this case satisfies Eq. (19) with the initial condition (24). 
As to the nuclear field $(j/, ri,r2), we may continue to consider it taken in the forward direction, 
neglecting the small nuclear momenta transferred to the gluons. So $ continues to depend only on 
ri2 and satisfies the same equation (26) and initial condition (42). 

Propagator P{y,ri,r2), apart from y depends on three variables for which we may take rf, 
and angle (f) between ri and r2- We introduce a grid in variables Inr^, Inr^ and <p, dividing the 
interval in the first two variables into A^i points and in (p into N2 points. The evolution equation (19) 
then requires calculation of the integrand at Nf ■ N2 points at each value of intermediate rapidities. 
This severely restricts the numbers of Ni and N2 admissible for given calculation facilities. As a 
result, in the past there were quite few calculations of the solution to the non-forward BK equation 
[19, 20] with modest values of Ni and 7V2. In our case we use the Runge-Kutta method of the 
solution with the maximal value Ni = N2 = 80 compatible with the reasonable calculation time. 
The limits in r were taken as before 10~^ < r < 10^. The initial condition (24) for the convoluted 
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Figure 14: Propagator P{y,ri = r2,(j) = '^/2) as a function of Y for the initial conditions (45) 
with ci = 10 and C2 = 2. Curves with maxima from right to left correspond to y = 2,4,6,8, 10 
respectively. 



propagator was taken at y' = in the form borrowed from [19] 



P{y 



(45) 



where 6 = (ri + r2)/2. We have studied 9 cases with c\ = 0.1, 1 and 10 and C2 = 0.2, 2 and 20. In 
all cases the results are quite similar. 

Calculations produce an array of data P{y, rf , r^, (p) on the three-dimensional grid. The results 
can be presented in various ways. In Fig. 14 we show values of P(?/, r^, r^, 0) with ci = 10 and 
C2 = 2 at ri = r2 and (p = Tr/2 for different y. 

As we see the general trend of the curves is the same as for the forward case. Fig 12, although 
the decrease of the height of the maximum starts from larger rapidities. The behaviour of P{y, ri = 
r2, (p) at different values of (f) and for different ci and C2 in the initial conditions (24) is the same, 
although the value of the rapidity Y at which the decrease starts is slightly different for different 
initial conditions. 

To have an overall picture of the behaviour of the non-forward propagator with rapidity, as for 
the forward case, we studied the convolution of three propagators with the triple pomeron vertex 
(dropping operator K, which hopefully does not change the behaviour at large y, see [8]) 



(fri(Pr2(Pr3 

''l2'"23''31 



Piiy, ri,r2)P2{y, r2,n)P2,{y, ra, ri). 



(46) 



where Pi, P2 and P3 may come from the same or different initial conditions. We again present 



/2(y) = a2e-^2(2/)y^ y = ay 



(47) 



and in Fig. 15 show values of A2(y) for four cases, with the three propagators taken 1) all for 
the initial condition (IC) with ci = 10, C2 = 2 (IC=A) 2)all for the IC with ci = 1, C2 = 20 
(IC=B) 3)two propagators with IC=A and one with IC=B and 4) two propagators with IC=B 
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Figure 15: Exponent A2 in (46) as a function of rapidity for the convolution of solutions to Eq. 
(19) with the initial conditions (IC) A and B in (45) (see explanation in the text. The curves from 
top to bottom correspond to 1) all three solutions with IC=A 2) two with IC=A and one with B 
3) two with IC= B and one with A and the lowermost to all three solutions with IC=B 

and one with IC=A. We see that at low values of Y the exponents A2 are different for these cases 
and predominantly negative, which means that the propagators in fact grow at these rapidities. 
However at larger rapidities different A2 converge to a common positive value indicating that the 
convolution decreases with rapidity and that this decrease is universal. 

6 Conclusions 

We have studied numerically the BFKL pomeron propagator in the external field created by the 
solution of the BK equation in the nuclear matter. We have found that for more or less arbitrary 
set of initial conditions the convoluted propagator vanishes at large rapidities, its coordinate de- 
pendence in the forward case tending to the (5-function. This gives reasons to believe that the 
propagator itself vanishes at large rapidities in the nuclear background. This result follows only 
with the field being the exact solution of the BK equation. 

Our results may be a starting point for the perturbative calculation of loop contributions in the 
nuclear background. However the technical difficulties seem to be quite formidable, since calcula- 
tions are only possible numerically and even simplest loops involve quite a number of spatial points 
over which one has to integrate. One needs to invent a reliable approximate representation for the 
non-forward propagator in the nuclear field which would allow one to perform all these integra- 
tions effectively. A somewhat simpler application is to the study of double inclusive cross-section 
for gluon jet production in nucleus-nucleus collisions, which only requires the forward propagator. 
Both these complicated problems are left for future investigations. 
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